Spectral functions of the Falicov-Kimball model with electronic ferroelectricity 
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We calculate the angular resolved photoemission spectrum of the Falicov-Kimball model with 
electronic ferroelectricity where d- and /-electrons have different hoppings. In mix- valence regimes, 
the presence of strong scattering processes between d-f excitons and a hole, created by emission 
of an electron, leads to the formation of pseudospin polarons and novel electronic structures with 
bandwidth scaling with that of d-f excitons. Especially, in the two-dimensional case, we find that 
flat regions exist near the bottom of the quasiparticle band in a wide range of the d- and /-level 
energy difference. 

PACS numbers: 71.27.-|-a,71.28.-|-d,77.80.-e 



Ferroelectric materials have long been important to 
physical research and technological application. Besides 
commonly known displacive and order-disorder mecha- 
nisms recently the idea of electronic ferroelectric- 
ity (EFE) has been attracting considerable attention. 
EFE was introduced by Portengen et al. 0, Q in their 
mean field theory of an extended Falicov-Kimball model 
(FKM) 0]. The FKM, having a long successful history in 
dealing with correlated electron systems, was then used 
to model a system consisting of two bands of different 
parity, say itinerant d electrons and localized / orbitals, 
in which mix-valence states may occur depending on the 
d-f coupling strength. Portengen et al. pointed out that 
inclusion of a d-f hybridization in the FKM could give 
rise to a spontaneous electric polarization due to a Bose- 
Einstein condensation (BEG) of d-f excitons when the 
excitation energy goes to zero at the critical value of the 
/-level energy y{. It is expected that such a purely 
electronic mechanism would provide fascinating physical 
features which are desirable for many applications; for 
instance, the static dielectric constant in an electronic 
ferroelectric could exceed 10'' 0. 

However, EFE in the FKM is still controversial after 
tested by different theoretical treatments of the under- 
lying strongly correlated electron system 0, 0, 0] . Most 
recently, Batista proposed a new extension of the FKM 
in which an /-/ hopping was included |^. Mapping the 
strong coupling limit of this model into an xxz pseu- 
dospin 1/2 model with a magnetic field along the z-axis 
and supported by quantum Monte Garlo calculations per- 
formed earlier by Schmid eJ al. , he showed that a BEG 
of d-f excitons does exist in the phase diagram 

The purpose of this Letter is to present a theoretical 
study of the angular resolved photoemission spectrum of 
this extended FKM. In the pseudospin picture 0, we 
study the spectral functions of one hole, created by emis- 
sion of one electron, in an extended t-J model with differ- 
ent hoppings for each pseudospin flavor. The problem of 
single hole motion in a local (pseudo)spin background has 



become an essential issue in understanding anomalous 
physical properties of high temperature superconductors 
;1(1 .Ih .12.1 and of colossal magnetoresistance mangan- 
ites [13, [l^. It has been demonstrated that quantum 
antiferromagnetic (pseudo)spin fluctuations have strong 
impact on the low-energy scale physics of these materials. 
In this Letter, we show that strong scattering processes 
between the hole and the d-f excitons take place even in 
the ferroelectric regime that has a ferromagnetic pseu- 
dospin configuration, thus give rise to the formation of 
pseudospin polarons and novel electronic structures. 

The extended Falicov-Kimball model for spinless 
fermions on a hypercubic lattice is 

i i (ij) (ij> i 

(1) 

where nf = d^dj and n( = fl f-^ are the occupation num- 
bers of d- and /-orbitals, respectively. Representing the 
two orbital flavors by a spin 1/2 variable, Ci| = di and 
Ci^ = /i, and using the notion of the pseudospin operator 
= E;.^ c\t,^i.uCiy with {ct^^} being the PauH matrices, 
one can derive an effective t-J model, iTeff = Ht -\- Hj, 
in the strong coupling limit of Q near half-filling 

(iJ> i 

Ht = +4^CiJ, (2) 

where q^. — Ci^(l — rii^) is the constrained fermion 
operator, t^ = ta, ti ^ tj , = 2{t] + tl)/Uf'^, 
J± ~ At-iti/U^'^, and = £d — £/ is the d- and /-level 
energy difference acing as a magnetic field along the z 
direction. 

At half filling only Hj contributes to Heg. The ground 
state of Hj is one of the following phases: full d-band, 
full /-band, staggered orbital ordering, and the BEG of 
d-f excitons that is ferroelectric for J± < or antiferro- 
electric for J± > We assume that the BEG of d-f 
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FIG. 1: Schematic pictures of the pseudospin configuration in 
(a) full /-band, (b) staggered orbital ordering, (c) ferroelectric 
canting state, and (d) antiferroelectric canting state. 
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FIG. 2: Two dimensional phase diagram of Hj obtained from 
linear spin wave theory (solid lines) and from quantum Monte 
Carlo calculations |3, j^l (dashed lines) . Shadow area denotes 
possible phase separation. 



excitons is single-phased - in the present case, the corre- 
sponding pseudospin configuration is a canting state with 
respect to the z axis [Fig.^c) or^fd)], where the canting 
angle 9 is determined by 



IB. 



zS{J, + \J^\)' 



(3) 



where S — 1/2 is the value of the pseudospins and z 
is the coordination number. Indeed, the spin excitation 
(i.e., d-f exciton) in the canting state is a gapless Gold- 
stone mode which corresponds to a uniform precession 
of the pseudospins around the z axis. Now we diagonal- 
ize Hj in linear spin wave theory and obtain the follow- 
ing phase diagram as shown in Fig. |5] The borderline 
between staggered orbital ordering and the BEG of d-f 
excitons is \Bz\ = zS^J J| — J\ and the corresponding 
phase transition is a first-order transition; the borderline 
between the BEG of d-f excitons and full d- or /- band 
is \Bz\ = zS{Jz + \J±\) and the corresponding transi- 
tion is a second-order transition. In a whole, our results 
agree with the quantum Monte Garlo studies on finite size 
systems up to 96 x 96 Notice that in Fig. |21the bor- 
derline between staggered orbital ordering and the BEG 
of d-f excitons obtained from linear spin wave theory 
(solid line) does not perfectly coincide with that from 
the Monte Garlo calculations (dashed line). This might 
result from our single phase approximation, indicating 
the existence of a phase separation regime between these 
two lines (shadow area). Nevertheless, the difference is 
small if the values of Jz and | J_l| are close to each other, 
and this is favorable for ferroelectricity to take place 
We shall analyze the hole dynamics in each phase in the 
following. 

Full f-band |l5| . The pseudospin configuration is a 
ferromagnetic state with all pseudospins down [Fig^Ia)]. 



The pseudospin wave spectrum is cj^"" ~ — zS{Jz -\- 

\J±\lci) where 7q = s'**^/-^ ^^'^ (5 is a unit vector 
connecting nearest neighbors. Since the hole motion will 
not change the pseudospin configuration, the hole can 
move freely with dispersion ej^"" = — zi|7k. 

Ferroelectric phase. The pseudospin configuration is 
considered here as a uniform canting state [Fig^c)]. The 
pseudospins are parallel, but they take new equilibrium 
direction which are tilted by a certain angle 9 to the z 
axis. This new axis of quantization is considered to lie 
in the x-z plane as a result of spontaneous symmetry 
breaking. We perform a uniform rotation of the orbitals 
about the y axis by 9, 



CiT = cos-ei| -sm-Cij^, Ci| = sm-ei| +cos-ei4, 
rf = S? cos 9 - S? sin 9, ^ S? cos 9 + S? sin 9, (4) 

to obtain the ferromagnetic configuration 
\ ■ ■ • S iS f^^- ■ ■) = ii • • •) as the vacuum state. 
Then, we employ the slave-fermion formalism to 
cop e with the constraint of no doubly occupancy 
[ll|. Defining holon (spinless fermion) operators hi 
so that = ^Jflj, — where Cj = S7 is the 
hard-core boson operator, we arrive at an effective 
pseudospin-polaron Hamiltonian in the momentum 
space 

Ht = J2^^^l^k + ^Mkcihlh^_qaq + H.C., (5) 

k k,q 

where a^s are pseudospin wave operators, flq = u^aq 4- 



UqaLq, with dispersion lu^ = : 
The transformation coefficients are u, 
Chl)-'^' + l]/2}^/' and vq = -sgn(C«7q){ 



q = {[-^e,q(^l,q ~ 
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Cg7q) — l]/2}^/^. Here the shorthand notations are 
Be = {-\J±\cosH + J.sinH - \J±\)/2, Ce = Be + \J±\, 
and Ig^q — \J±\ + Bejq- Hence ujq=o = 0. The bare 
hole dispersion is Ej^. = —z{ti cos^ 5+^1 ^i'^^ i)7k, and 
the hole-pseudospin-wave couphng is 



Ml 



kq 



(7k-qWq +7kWq)- (6) 



When both pseudospin flavors have the same hopping in- 
tegral, i.e., i| = = t, Ht has the SU(2) symmetry, lead- 
ing to a vanishing Mkq and a free hole propagation with 
dispersion — zt7k. Hence, the fact that t-\ ^ ti accounts 
for the formation of the pseudospin-polaron even in the 
presence of a ferromagnetic pseudospin background. 

Using the self-consistent Born approximation (SCBA) 
in which the spectral functions of one hole in a t-J-like 
model can be accurately calculated jO, 0, lisl , we 
compute the hole Green's function G(k, w) = [u — e^^ — 
S(k, a;) + iO+]~^ self-consistently with the self-energy. 



E(k,c.) 



E^^kqG(k 



q,w - 



(7) 



Thus, the spectral functions are given by A(k, uj) — 
— ImG(k, CL')/7r, and the quasiparticle dispersion is i?k = 

ek + ReS(k,£;k). 

Antiferroelectric phase. This phase [Fig^d)] is possi- 
ble when Jj_ > 0. We find that linear spin wave theories 
of the antiferroelectric and ferroelectric canting states are 
formally connected with a displacement of momentum: 
q — !■ q — Q where Q —{n, . . . ,Tr). The effective Hamilto- 
nian for hole hopping is 



Ht^Yl ^k^lc'ik + ^^kq/li^lk-qaq-Q + H.< 

k k,q 

where = —z{ti cos^ f ~ ^i'^^ f )7k and 



(8) 



Afkq = — sin t 



(7k-qUq-Q + Tkl'q-q). (9) 



Likewise, we employ (fTj) to calculate the hole self energy. 

Staggered orbital ordering. The pseudospin configura- 
tion is a z-directional Neel state [Fig 1(b)]. Considering 
the Neel state as the vacuum state, we define holon op- 



erators hi and gj so that 



h 



/itfl; on the I 

sublattice and Tiji = gjb^, Cj| = gj on the t sublattice. 
Here li — t7 on the J, sublattice and 6j = on the 
t sublattice are hard-core boson operators. We arrive at 
an effective pseudospin-polaron Hamiltonian 

Ht = E 'Afkq/lk5k-qaq + ^kqSk'ik-q/^q + 
k,q 

where aq and /3q arc pseudospin wave operators, 

flq = Uqaq -|- VqP'Lq and 6_q = UqP-q + 
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FIG. 3: (a) presents the quasiparticle dispersion (solid 
line, left scale), and pseudospin wave dispersion Uq (dashed 
line, right scale) in the ferroelectric phase with Bz = 0.2. (b)- 
(d) shows the hole spectral functions A(k,u}) along (b) (0, 0) — 
(7r,7r), (c) (tTjTt) — (tt, 0) and (d) (tt, 0) — (0,0) directions. 



Bz, respectively, and Uq — {[zSJz/{i.^q + ui^) + 

l]/2}l/2, ^ -Sgn(Ji7q){[^^J./« + ^q) " 

l]/2}^/^. The hole-pseudospin-wave coupling is 

Mkq = -2;^2/A^(t|7k-qUq -|- tlJkVq) and Lkq = 

— z y/2/ N {t ijk-qUq + t|7kWq). Thc summatious over k 
and q are restricted inside the reciprocal Brillouin zone 
of one sublattice. Within the SCBA, we obtain the fol- 
lowing self-consistent equations for the two types of hole 
propagators: Gj(k, w) = [uj — Ej(k, w) -f- iO+J"-^ where 
j = 



WqaJ- , with dispersion 



zS{J! - Ji7^)V2 ± 



= g,h with E/,(k,a;) = J2q' M^qGgCk- (i,uj ~ uj+) and 
(k,w) = Y.q'LlqGhO<i~ ci, ui-u)-). 
To facilitate EFE, the system must be in a mix-valence 
regime and the two bands involved should have differ- 
ent parity, thus t| and t^ have opposite signs, so do 
Jz and J±; furthermore, it would be more probable 
if both bands have similar bandwidths 0. Hence, we 
adopt the following parameters 0|: i| = 1, = —0.8, 
U^''^ = 10, hence Jz = 2(f2 + i2)/j//rf ^ o.328 and 

= At^tjU^'^ = -0.32, with Bz = Cd - Cf being a 
free parameter ~- in a real material, changing B^ could be 
achieved by applying pressure or alloying. Our numer- 
ical results are illustrated in the two-dimensional case. 
In this case, the system transits from staggered orbital 
ordering to the BEC of d-f excitons at \Bz\ — 0.144 and 
then to a non- mixed- valence regime at \Bz\ — 1.296 as 
Bz increases. All numerical calculations were carried out 
on a 16 X 16 square lattice. 

In Figs. 01 El andlHl we show the spectral functions and 
the quasiparticle dispersions in the ferroelectric regime 
with Bz = 0.2, 0.8, and 1.2, respectively First of aU, 
as presented in Figs. Ol^a), 2{a), and^a), the quasipar- 
ticle bandwidth W and the pseudospin wave bandwidth 
I^spw have similar values. Such band narrowing can be 
understood in the following way: gapless pseudospin ex- 
citations are easily stimulated by incoherent hole motion. 
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FIG. 4: The same quantities as in Fig. |3]but with = 0.8. 

leading to the formation of the quasiparticle (QP), pseu- 
dospin polaron, which is the propagating hole surrounded 
by a cloud of polarized pseudospin waves. Therefore, W 
does not scale with hopping integrals but with IVspw 

As displayed in Figs.|3Ib)-(d), for any k, there is a well- 
defined quasiparticle pole (i.e., zero pseudospin wave) 
at the low energy side which is well separated from a 
broad, incoherent, multiple-pseudospin-wave background 
extending to the full free-electron bandwidth. Fig. |3ta) 
shows that the bottom and the top of the QP band locate 
at (7r/2,7r/2) and (0,0), respectively. Furthermore, there 
is a flat region around (tt, 0) at which the QP energy 
is close to its minimum, leading to a strongly distorted 
density of states with a massive peak near the bottom 
of the QP band. These features obtained for single hole 
motion in the ferromagnetic pseudospin background are 
quite similar to those in the cuprate t-J model where the 
hopping integrals for both spin flavors are the same and 
the spin background is a Heisenberg antiferromagnet |l2l | . 
Note that the pseudospin wave spectrum does not vanish 
at (tt, tt) while it is gapless at (0,0). Consequently, the 
QP energy at (tt, tt) is lower than that at (0,0). 

Further increasing will continue lowering the QP 
energy at (7r,7r). When = 0.8, (7r,7r) becomes the 
QP band bottom, as shown in Fig. 2{a) . While the flat 
region around (tt, 0) is shrinking, a new, surprising flat 
region appears around (tt, tt). Therefore, the feature of 
an extended van Hove singularity near the bottom of the 
QP band survives up to B^ = 0.8. In Figs. IHb)-(d), the 
spectral functions are still characterized by sharp peaks 
at the low energy side which is well separated from a 
broad, incoherent background. However, apart from the 
QP band bottom, considerable spectral weights move to 
the incoherent part, indicating that the hole motion is 
more severely damped by incoherent processes. 

When Bz is increased to 1.2 [Figs. |3b)-(d)], almost 
all spectral weights locate in the incoherent part ex- 
cept near the QP band bottom, (tt, tt). This means that 
the hole spectral weights ZCk) = [1 — dT,{k,uj)/duj]~^^^ 
are strongly reduced by a cloud of polarized pseudospin 
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FIG. 5: The same quantities as in Fig.|3]but with = 1.2. 

waves. The shape of the QP dispersion [Fig. [S^a)] be- 
comes similar to that of free hole dispersion (x 7k, yet 
with bandwidth W ~ Wgpw As B^ further increases, 
Z(k) at wave vectors not very close to (tt, tt) goes to zero 
quickly until the pseudospin polaron picture becomes in- 
valid at Bz = 1.296. Then, the system evolves into full 
/-band with free hole motion, implying that ZCk) = 1. 

To summarize, we present a systematic study on the 
evolution of the spectral functions of a single hole in 
the FKM as a function of B^, the d- and /-level energy 
difference. We find that strong scattering processes be- 
tween the hole and d-f excitons exist in the mix-valence 
regimes, thus lead to the formation of pseudospin po- 
larons and interesting electronic structures with band- 
width scaling with that of the pseudospin excitation, even 
for single hole propagation in a ferromagnetic pseudospin 
background (i.e., the ferroelectric phase) due to the dif- 
ferent hopping integrals of the d- and /-electrons. Fur- 
thermore, flat regions near the bottom of the quasipar- 
ticle band are found in a wide range of B^ in the two- 
dimensional case. As we know, similar anomalous, yet 
crucial hole spectral features were found in high temper- 
ature superconductors 0, 0|, thus we speculate that 
besides EFE, the extended Falicov-Kimball model could 
have many novel properties, especially upon hole doping. 
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